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Among  other  applications,  the  study  of  thermal  properties  of  large  networks  of  carbon  nanoparticles 
may  have  a  critical  impact  in  loss-free ,  more  compact  and  efficient  thermal  storage  systems,  as  well  as 
thermally  conducting  polymeric  materials  for  innovative  low-cost  heat  exchangers.  In  this  respect,  here, 
we  both  review  and  numerically  investigate  the  impact  that  nanotechnology  (and  in  particular  carbon- 
based  nanostructures)  may  have  in  the  near  future.  In  particular,  we  focus  on  the  role  played  by  some 
geometrical  and  chemical  parameters  on  the  overall  thermal  transmittance  of  large  complex  networks 
made  up  of  carbon  nanotubes  (CNTs),  that  can  be  potentially  added  as  fillers  to  (thermally)  low- 
conductive  materials  for  enhancing  the  transport  properties.  Several  configurations  consisting  of  sole 
and  pairs  of  single-walled  carbon  nanotubes  (SWCNTs)  and  double-walled  carbon  nanotubes  (DWCNTs), 
characterized  by  different  dimensions  and  number  of  C-O-C  interlinks,  are  considered.  Based  on  the 
results  found  in  the  literature  and  using  focused  simulations  using  standard  approaches  in  Non- 
Equilibrium  Molecular  Dynamics  (NEMD),  we  highlight  the  dependence  on  the  particle  diameter,  length, 
overlap  and  functionalizations  of  both  thermal  conductivity  and  boundary  resistance  across  CNTs,  which 
are  indeed  the  relevant  quantities  for  obtaining  composite  materials  with  desired  unusual  thermal 
properties.  We  observe  that  CNTs  with  short  overlap  length  and  a  few  interlinks  already  show  a 
remarkable  enhancement  in  the  overall  transmittance,  whereas  further  increase  in  the  number  of  C-O-C 
connections  only  carries  marginal  benefits.  We  believe  that  much  understanding  has  been  gained  so  far 
in  this  field  thanks  to  the  work  of  chemists  and  material  scientists,  thus  it  is  time  to  draw  the  attention  of 
engineers  active  in  the  energy  sector  and  thermal  scientists  on  such  findings.  Our  effort,  therefore,  is  to 
gather  in  this  article  some  guidelines  towards  innovative  thermal  systems  that  may  be  manufactured 
and  employed  in  the  near  future  for  addressing  a  great  challenge  of  our  society:  Storage  and  use  of 
thermal  energy. 
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1.  Introduction 

Thermal  energy  storage  for  sufficiently  long  periods  is  known  to 
be  a  critical  issue  in  thermal  sciences.  In  general,  nanotechnology 
(and  particularly  carbon  based  nanostructures)  may  hold  the  premise 
for  successfully  addressing  energy  issues  in  modern  society  [1,2]. 

In  particular,  recently,  an  innovative  system  with  the  capacity  of 
storing  thermal  energy  potentially  for  an  indefinite  time  have 
attracted  much  attention  thanks  to  the  pioneering  work  of  scien¬ 
tists  at  the  Fraunhofer  Institute  for  Interfacial  Engineering  and 
Biotechnology,  together  with  ZeoSys  GmbH  [3  .  In  addition,  such  a 
new  technology  can  store  three  to  four  times  more  energy 
compared  to  water  based  systems.  The  new  system  contains  zeolite 
pellets  and  one  of  the  key  challenge  here  is  (i)  to  ensure  optimal 
heat  transfer  (high  overall  conductivity)  through  the  zeolites 
during  the  thermal  loading  and  (ii)  to  allow  water  percolation 
during  the  thermal  release,  in  the  same  device.  The  latter  compet¬ 
ing  (and  extreme)  features  may  be  met  by  dispersing  highly 
conductive  materials  (fillers)  within  the  zeolite  pellets  (matrix)  [4  . 

The  use  of  carbon-based  nanoparticles  (e.g.  carbon  nanotubes, 
graphene,  etc.)  certainly  represents  one  of  the  most  promising 
alternative.  On  the  other  hand,  even  for  perfectly  dispersed  fillers 
forming  an  ideal  percolating  network,  realizing  a  sufficiently  high 
thermal  transmittance  for  the  overall  network  is  everything  but  a 
simple  task.  In  fact,  typically  high  thermal  resistances  are  experi¬ 
enced  by  the  heat  flux  across  the  several  interfaces  characterizing 
such  composite  materials  (see  e.g.  [5,6  ).  It  is  therefore  highly 
desirable  to  systematically  investigate  on  the  several  mechanisms 
that  contribute  to  the  value  of  the  overall  thermal  transmittance  of 
such  complex  percolating  networks.  Moreover,  those  studies  have 
a  critical  impact  also  on  polymeric  materials  for  innovative  heat 
exchangers  and  other  composite  materials  [7-10  . 

Since  their  first  observation  by  Iijima  [11],  Carbon  Nanotubes 
(CNTs)  have  attracted  an  increasing  attention  because  of  their 
outstanding  physical  properties.  In  fact,  superior  mechanical, 
electric  and  (most  importantly)  thermal  properties  of  CNTs  have 
been  reported  [12],  and  this  certainly  makes  such  material 
a  versatile  building  block  for  future  engineering  applications  [13]. 

Among  others,  thermal  properties  of  CNTs  have  been  widely 
investigated  in  recent  years,  both  experimentally  and  theoretically 
[14,15  .  For  instance,  pristine  CNTs  are  characterized  by  one  of  the 
largest  value  of  thermal  conductivity,  ranging  from  2400  W  m”  K“ 
to  4740  Wm-1  K“  depending  on  the  experimental  conditions  [14  . 
Moreover,  thermal  conductivity  of  CNTs  may  be  tuned  by  introducing 
defects,  vacancies,  surface  functionalizations  or  junctions  [16,17],  thus 
allowing  the  design  of  CNT-based  composite  materials  with  desired 
enhanced  thermal  properties  [18,19  .  CNT-based  composites  are  typi¬ 
cally  made  out  of  polymeric,  ceramic  or  metallic  matrices  and  carbon- 
based  fillers,  such  as  CNTs,  graphene  or  carbon  nanoribbons  [20,21]. 

CNTs  spontaneously  tend  to  bundle  and  cross-link  each  other 
because  of  their  high  intermolecular  cohesive  forces  [22  ,  whereas 
precise  geometries  of  CNT  networks  can  be  induced  by  means  of 
proper  chemical  functionalizations  [23  .  Therefore,  CNT  networks 
are  energetically  and  thermodynamically  stable  up  to  600  K,  and 
CNTs  synthetized  within  a  polymeric  or  ceramic  matrix  can  exist  in 
a  stable  and  designable  cross-linked  form  [24  .  Due  to  their 
peculiar  thermal  and  technological  properties,  CNT  fillers  have 
been  employed  in  a  promising  variety  of  applications,  such  as 
thermal  interface  material  [25-27],  electrical  interconnects  [28], 


electronic  nanocircuits  [29],  triggered  drug  release  [30],  artificial 
muscles  [31],  thermal  nanocircuits  [32-34],  nanofluids  for  heat 
transfer  [35-37],  energy  harvesting  [38]  and  gas  or  thermal 
storages  [24,39].  For  example,  CNT  networks  are  promising  tech¬ 
nological  solutions  for  improving  the  internal  distribution  of  heat 
within  thermal  storage  applications,  such  as  nanoporous  zeolite 
devices,  in  order  to  reach  a  better  control  and  efficiency  of  the  heat 
accumulation/release  cycle  [40  .  In  the  latter  case,  CNT  networks 
need  to  be  engineered  and  properly  designed  for  reducing  thermal 
losses  at  low  filling  ratio  in  the  matrix.  However,  experimental 
benchmarks  and  molecular  simulations  of  carbon-based  networks 
reveal  an  overall  thermal  performance  lower  than  expected,  which 
is  mainly  due  to  CNT-CNT  and  CNT-matrix  thermal  boundary 
resistances  [41-44]. 

It  is  well  known  that  heat  conduction  by  CNTs  is  predominantly 
a  phononic  phenomenon  (i.e.  coupled  vibration  of  carbon  atoms), 
rather  than  an  electronic  one  [45].  The  reasons  for  high  thermal 
conductivity  of  CNTs  are  due  to  their  large  speed  of  sound,  long 
range  crystallinity  and  long  phonon  mean  free  path  [14,46].  Never¬ 
theless,  the  bottleneck  of  heat  transfer  within  carbon-based  net¬ 
works  is  the  thermal  boundary  resistance  (also  known  as  Kapitza 
resistance),  which  defines  the  thermal  impedance  between  differ¬ 
ent  materials.  Kapitza  resistance  is  generated  by  phonon  scattering, 
which  produces  a  local  temperature  discontinuity  in  heat  transfer 
phenomena  at  the  nanoscale.  Vacancies,  functionalizations,  defects 
as  well  as  interfaces  between  different  materials  (solid  or  liquid) 
contribute  to  thermal  boundary  resistance  [47-55  .  Therefore,  the 
difference  between  intertube  and  intratube  thermal  conduction, 
which  can  also  be  orders  of  magnitude,  limits  the  overall  thermal 
transmittance  of  the  CNT-based  composites.  Hence,  a  better  under¬ 
standing  of  the  geometry,  chemistry  and  arrangement  of  CNT 
networks  is  needed,  in  order  to  have  guidelines  in  designing 
composite  materials  with  desired  heat  transfer  performances. 
Recent  studies  about  carbon-based  composite  materials  show  that 
a  better  efficiency  in  heat  transfer  is  achieved  when  CNTs  percolate 
within  the  matrix,  creating  a  homogeneous  and  cross-linked  net¬ 
work  [6,56,57  .  Owing  to  the  current  experimental  difficulties  in 
synthetizing  high  quality  and  well-ordered  CNTs,  computational 
studies  have  been  conducted  to  systematically  analyze  the  thermal 
properties  of  CNT  and  CNT  networks  [57-60].  Among  others, 
Molecular  Dynamics  (MD)  simulations  have  been  widely  used  for 
improving  the  physical  understanding  and  predicting  the  thermal 
transport  in  carbon-based  composite  materials  [61  .In  this  study, 
thermal  properties  of  CNTs  with  different  diameter,  length,  chir¬ 
ality,  number  of  walls,  overlap  and  chemical  functionalization  are 
investigated  by  means  of  MD  simulations.  Hence,  this  work  aims  at 
systematically  reporting  and  analyzing  the  geometrical  and  che¬ 
mical  parameters  affecting  the  thermal  conductivity  and  the 
thermal  boundary  resistance  of  CNT  networks,  in  order  to  come 
up  with  a  few  general  guidelines  for  designing  the  overall  thermal 
transmittance  of  CNT-based  composites. 


2.  Materials  and  methods 

Non-Equilibrium  Molecular  Dynamics  (NEMD)  simulations  are 
performed  for  estimating  the  thermal  conductivity  X  of  CNTs  and 
the  thermal  boundary  resistance  Rk  between  adjacent  CNTs.  The 
latter  configuration  can  be  regarded  as  the  building  block  of  large 
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and  complex  carbon-based  networks  as  those  of  interest  in  today 
nanofluids  and  composite  material  [5,7,62]. 

Thermal  properties  of  Single-Walled  Carbon  Nanotubes  (SWCNTs, 
Fig.  la)  and  Double-Walled  Carbon  Nanotubes  (DWCNTs,  Fig.  lb), 
which  are  composed  by  two  coaxial  SWCNTs,  are  investigated  in 
different  setups,  ranging  from  a  sole  CNT  to  larger  CNT  aggregates.  First 
of  all,  the  effect  of  the  CNT  geometry  on  thermal  conductivity  is 
investigated.  In  particular,  in  our  simulations  we  choose  (i)  CNT  with 
diameters  of  0.7, 1.4,  4.5  and  5.0  nm;  (ii)  CNT  with  length  of  10,  20,  30 
and  40  nm;  (iii)  CNT  with  different  chirality,  namely  armchair  (m,m), 
chiral  (m,n)  or  zigzag  (m,0).  Second,  the  effect  of  tube  overlap  (ranging 
from  0%  to  100%  of  the  CNT  length)  and  chemical  bonding  (0-8  C-O-C 
joints  between  parallel  CNTs,  Fig.  lc  and  d)  on  thermal  boundary 
resistance  between  contiguous  CNTs  is  explored. 

In  our  simulations,  the  reciprocal  orientation  between  the  CNTs 
jointed  by  oxygen  atoms  is  initialized  according  to  C-O-C  bond 
structure  (Fig.  le,  C-0  bonds  are  1.43  A  long,  C-O-C  angle  109.5° 
wide)  [63],  while  two  types  of  interactions  are  taken  into  account: 
(i)  the  bonded  interactions  modeled  as  harmonic  stretching  and 
angle  potentials,  whereas  the  Rickaert-Bellemans  potential  is  used 
for  proper  dihedral  angles;  (ii)  non-bonded  interactions  accounting 
for  the  Van  der  Waals  forces  and  modeled  as  a  12-6  Lennard -Jones 
potential  with  a  shifted  cut-off  function  (1.2  nm  cut-off  distance). 
All  details  on  the  employed  force-field  parameters  are  reported 
elsewhere  [5,62-65]. 

The  equations  of  motion  are  solved  by  a  leap-frog  algorithm 
with  time  step  of  At=  0.5  fs;  whereas  periodic  boundary  condi¬ 
tions  along  x,y,z  directions  are  imposed.  Upon  energy  minimiza¬ 
tion,  the  system  is  initialized  at  300  K  (Maxwellian  distribution  of 
velocities)  and  coupled  to  a  Berendsen  thermostat  for  50  ps,  until 
system's  energies  relax  to  a  steady  state.  Afterwards,  the  NEMD 
approach  is  used  for  evaluating  the  thermal  properties  of  CNT  or 
CNT  network  [66  .  A  temperature  gradient  is  induced  in  the  CNT 
(or  in  the  CNT  network)  by  means  of  two  Nose-Hoover  thermo¬ 
stats  [67,68]  (Ti  =  280  K;  T2  =  320  K)  attached  to  the  carbon  atom 
rings  at  their  ends  (Figs.  la-d).  Thermostats  are  introduced  in 
different  positions  according  to  the  thermal  properties  to  be 
evaluated.  For  example,  in  the  DWCNT  setup  in  Fig.  lb  the  green 
dotted  thermostated  carbon  rings  are  employed  for  evaluating  X  of 


the  outer  CNTs,  whereas  the  red  ones  are  used  for  X  of  the  inner 
one.  Instead,  in  the  case  of  overlapping  SWCNTs  depicted  in  Fig.  lc 
the  green  dotted  thermostated  carbon  rings  are  employed  for 
evaluating  X  of  the  lower  SWCNT,  whereas  the  red  ones  for 
calculating  the  interlayer  Kapitza  resistance  between  the  SWCNTs. 
The  amount  of  atoms  in  the  canonical  ensemble  (NVT)  is  optimized 
for  obtaining  constant  heat  flux  and  convergent  values  of  X  in  the 
simulated  setups,  and  it  is  properly  scaled  according  to  the  system 
size  [69];  the  remaining  atoms  are  in  the  microcanonical  ensemble 
(NVE).  After  the  energy  minimization,  the  first  and  the  last  rings  of 
carbon  atoms  are  restrained,  in  order  to  have  fixed  boundary 
conditions  (Figs,  la-d)  [70  .  NEMD  simulations  are  continued  up 
to  4  ns  (depending  on  the  system  size).  Temperature  profiles  and 
(Nose-Hoover)  friction  factors  [5,67-69]  are  then  extracted  from 
the  simulated  trajectories,  in  order  to  evaluate  temperature  gra¬ 
dient  and  specific  heat  flux  in  the  NVE  region,  respectively  [5  . 
Steady  state  is  reached  when  the  thermal  conductivity  X  and 
boundary  resistance  Rk  (evaluated  every  200  ps)  tend  to  an 
asymptotic  value;  this  is  typically  achieved  after  half  of  the 
simulation  time  has  elapsed.  Finally,  X  and  Rk  values  are  evaluated 
over  a  data  range  of  at  least  500  ps  once  steady  state  is  achieved, 
and  checked  for  energy  fluctuations. 

All  simulations  are  performed  using  the  software  package 
GROMACS  [71],  while  picture  rendering  is  accomplished  using 
UCSF  Chimera  [72]. 


3.  Theoretical  background:  Thermal  transmittance  of  carbon 
nanotube  networks 

Thermal  conduction  in  CNT-based  composites  (e.g.  carbon 
fillers  dispersed  within  a  polymeric  and/or  ceramic  matrix)  is 
mainly  affected  by  the  mechanism  of  heat  diffusion  (i)  through 
CNTs,  (ii)  between  the  CNTs  in  a  CNT  network,  (iii)  between  the 
CNT  network  and  the  matrix,  and  (iv)  through  the  matrix.  Heat 
insulation  or  heat  conduction  enhancement  in  carbon-based  com¬ 
posite  materials  can  be  influenced  by  varying  the  amount  (and 
type)  of  cross-interactions  as  well  as  the  arrangement  of  the  CNTs 
within  a  network,  such  as  choosing  between  randomly  oriented 


Fig.  1.  Setups  considered  in  the  computations:  (a)  single-walled  carbon  nanotube;  (b)  double-walled  carbon  nanotube;  (c)  a  pair  of  overlapping  SWCNTs  connected  by  4 
C-O-C  covalent  bonds  (links);  (d)  a  pair  of  overlapping  DWCNTs  connected  by  4  C-O-C  covalent  bonds  (links);  and  (e)  detail  of  C-O-C  joints.  (For  interpretation  of  the 
references  to  color  in  this  figure,  the  reader  is  referred  to  the  web  version  of  this  article.) 


M.  Fasano  et  al.  /  Renewable  and  Sustainable  Energy  Reviews  41  (2015)  1028-1036 


1031 


CNTs  or  aligned  ones  [14  .  The  aim  of  the  present  study  is  to 
systematically  investigate  the  more  fundamental  aspects  (i.e.  at 
those  acting  at  the  nanoscale  level)  involved  in  the  conductive  heat 
transfer  for  the  above  setup. 

Following  a  standard  approach  in  the  macroscopic  heat  transfer 
phenomena,  the  specific  heat  flux  cp  generated  by  a  temperature 
difference  AT  within  a  CNT  network  can  be  defined  as 


(p  =  u  AT, 


where  the  specific  thermal  transmittance  u  can  be  evaluated  as 


1 

u 


2  y+ +  2  Rk,n, 
i  Ai  ij  n 


4.  Results 


4.1.  Thermal  conductivity  of  carbon  nanotubes 


A  first  set  of  simulations  has  been  performed  for  systematically 
assessing  how  geometrical,  chemical  and  physical  factors  affect 
thermal  conductivity  of  SWCNTs.  Moreover,  the  adopted  MD  setup 
and  force-field  parameters  have  been  validated  by  comparing  the 
obtained  trends  for  thermal  conductivity  with  experimental  and 
computational  works  published  in  the  literature. 

According  to  the  NEMD  approach,  thermal  conductivity  is 
measured  by  means  of  the  Fourier's  law: 


(p 

dT / dx ’ 


with  Ai  being  the  thermal  conductivity  of  the  ith  CNT  (with  length 
Li);  Rkij  the  thermal  boundary  resistance  between  ith  and  jth  CNTs; 
Rk  n  the  nth  thermal  boundary  resistance  between  the  CNT  network 
and  a  matrix  [73,74],  which  is  here  considered  just  in  series  for  the 
sake  of  simplicity  (not  investigated  in  this  work). 

More  specifically,  the  present  MD  simulations  aim  at  evaluating 
and  Rk  ij,  which  are  critical  for  the  accurate  estimate  of  the  overall 
thermal  transmittance  in  CNT  networks  of  increasing  complexity.  Here, 
we  focus  on  (i)  sole  SWCNT  (Fig.  2a),  where  \/u  =  L/A,  L  is  the 
nanotube  length  and  A  its  thermal  conductivity;  (ii)  sole  DWCNT 
(Fig.  2b),  where  1  /u^L((l  /2in)  +  (l  /20Ut))+fi/<  in  the  specific  setup  here 
considered  (heat  flows  from  the  leftmost  end  of  the  inner  CNT  to  the 
rightmost  end  of  the  outer  CNT),  Ain  and  Aout  are  thermal  conductivities 
of  the  inner  and  outer  CNTs,  respectively,  and  Rk  is  the  Kapitza 
resistance  between  the  layers;  (iii)  two  adjacent  SWCNTs  (Fig.  2c), 
where  1  /u  =  (L\/Ai)+Rk,u  +  (L2/A2)  and  Rk  u  is  the  Kapitza  resistance 
between  the  SWCNTs;  (iv)  two  adjacent  DWCNTs  (Fig.  2d),  where 

1  /  U=  [Li  ((1  /M,in)  +  0  / A\  ,out))  +  Rk,  1  ]  +  Rk,U  +  [^2(0  /  ^2, in)  +  ( V  ^2,  out))  +  Rk,2] 

in  the  analyzed  configuration  (heat  flows  from  the  leftmost  inner  CNT 
to  the  rightmost  inner  CNT).  In  the  latter  case,  Ayin  and  Ay0ut  are  the 
inner  and  outer  thermal  conductivities  of  the  top  DWCNT;  A2jn  and 
A2lout  are  the  inner  and  outer  thermal  conductivities  of  the  lower 
DWCNTs;  Rk/l  and  Rk  2  are  defined  as  the  intralayer  Kapitza  resistances 
in  the  DWCNTs  and  Rk  u  is  the  interlayer  Kapitza  resistance  between 
the  DWCNTs.  Clearly,  different  heat  paths  within  the  network  may 
imply  different  estimates  of  1/u.  The  above  formulas  are  reported  here 
only  with  the  purpose  to  better  highlight  the  dependence  of  u  on  the 
most  relevant  quantities  to  be  investigated  below. 


where  cp  =  @/S ,  with  &  being  the  steady  state  heat  flux  between 
two  thermostated  regions  flowing  across  the  tube  cross-section  S. 
Starting  from  the  initial  state,  where  the  carbon  rings  at  CNT  ends 
are  thermostated  at  320  and  280  K  respectively,  while  the  remain¬ 
ing  rings  are  equilibrated  at  300  K,  the  temperatures  of  internal 
carbon  rings  tend  to  a  steady  state  in  a  few  ps  (Fig.  3).  The  slope  of 
time  averaged  temperature  profile  dT/dx  is  measured  along  the 
x-axis  of  the  nanotube  in  the  region  not  attached  to  thermostats 
(Fig.  3).  Different  ways  of  computing  S  can  be  found  in  the 
literature,  which  significantly  affect  the  absolute  values  of  A.  Here, 
the  cross-sectional  area  is  intended  as 

S  =  2  nrb,  (4) 


Fig.  3.  NEMD  approach  for  evaluating  thermal  conductivity  A  of  CNTs  by  means  of 
Fourier's  Law. 


Fig.  2.  Quantities  appearing  in  the  calculation  of  thermal  transmittance  in  the  studied  setups:  (a)  SWCNT,  (b)  DWCNT,  (c)  a  pair  of  overlapping  SWCNTs  connected  by 
4  C-O-C  covalent  bonds,  and  (d)  a  pair  of  overlapping  DWCNTs  connected  by  4  C-O-C  covalent  bonds. 
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where  r  is  the  CNT  radius  where  thermostat  (i.e.  the  heat  flux) 
is  applied  and  b  =  0.34  nm  is  the  Van  der  Waals  thickness  of  carbon 
atoms  [75]. 

(10,10)  SWCNT  and  (5,5)  SWCNT,  with  length  of  10  nm,  are 
considered  as  benchmarks  for  validating  the  present  MD  results 
(Fig.  la):  The  measured  thermal  conductivities  (78  Wm-1  K_1  and 
75  W  m“  K_1,  respectively)  are  in  good  agreement  with  other  results 
found  in  the  literature  for  same  setups.  For  example,  Alaghemandi 
et  al.  measured  79  and  69  W  m”  K-1,  respectively  [32].  Moreover, 
we  found  that  restraints  on  carbon  rings  at  the  SWCNT  ends  (blue  area 
in  Fig.  la)  have  negligible  effects  on  A,  being  76  W  m“  I<“  the 
measured  thermal  conductivity  of  (5,5)  SWCNT  10  nm  long  with  free 
ends.  It  is  well  known  that  thermal  conductivity  of  CNTs  computed  by 
the  MD  simulations  can  be  orders  of  magnitude  lower  than  the 
experimental  one,  because  of  the  strong  size  dependence  of  A.  The 
reason  for  that  is  the  phonon  scattering  at  boundaries  [32,76],  which 
occurs  when  CNT  length  is  smaller  than  the  phonon  mean  path  (about 
1.5  pm  [77]).  In  this  work,  though,  we  do  not  focus  on  the  absolute 
values  of  thermal  conductivity  of  CNTs,  rather  we  intend  to  investigate 
the  impact  of  geometry  and  chemical  functionalizations  (i.e.  C-O-C 
tube  interlinks)  on  A. 

Thermal  conductivity  of  sole  CNTs  is  affected  by  geometry  (i.e. 
diameter,  length,  chirality  and  amount  of  layers)  and  chemistry  (i.e. 
defects,  vacations  and  functionalizations),  as  well  as  by  their 
physical  conditions  (e.g.  temperature).  While  the  latter  physical 
conditions  are  illustrated  in  other  works  [76],  geometrical  and 
chemical  factors  are  of  main  interest  here. 

First  of  all,  the  effect  of  the  tube  diameter  on  A  is  evaluated  for 
10  nm  long  SWCNTs,  by  choosing  values  ranging  from  0.7  up  to 
5.0  nm.  Fig.  4a  shows  that  thermal  conductivity  scales  linearly  with 
the  diameter  D  in  the  considered  interval  (fitting  parameters: 
k\=  3.0  W  m_1  K_1  nm-1  and  k2  =  73.5  Wm_1  K_1;  coefficient  of 
determination:  R2  =  0.99),  namely  from  75  W  m-1  K_1  (D= 0.7  nm) 
to  88Wm-1K_1  (D=5.0nm).  Despite  the  fact  that  there  is  no 
general  agreement  in  the  literature  on  the  effect  of  CNT  diameter  on  A 
[14],  the  relationship  A  ~  D  is  for  example  also  found  by  Qiu  et  al.  [78  . 
The  concurrent  decrease  of  phonon  scattering  and  increase  of 
conduction  channels  with  tube  diameter  may  provide  an  explanation 
of  the  above  slight  direct  proportionality.  However,  if  the  CNT  cross- 
sectional  area  is  evaluated  as  S  =  nr2  instead  of  S  =  2jirb ,  the  D  —  A 
relation  changes  as  A  ~  1/D;  whereas  thermal  transmittance 
U  ~  IS  increases  with  CNT  diameter  in  both  cases. 


Moreover,  the  effect  of  CNT  length  is  assessed  by  considering  (5,5) 
SWCNTs  10-40  nm  long.  Results  range  from  75  (L=10nm)  to 
186Wm_1K_1  (L=40nm).  Trend  in  Fig.  4b  suggests  that  A~La 
(fitting  parameters:  /<3  =  19.2  W  m“  1  K"  nm"01  and  a=0.62;  coeffi¬ 
cient  of  determination:  R 2  =  0.99),  where  a  is  in  good  agreement 
with  other  works  [32].  Therefore,  the  increase  in  the  tube  length  has  a 
significant  effect  on  A:  In  the  ballistic  regime,  longer  CNTs  also  present 
vibrational  modes  with  longer  wavelengths,  which  add  new  path¬ 
ways  for  phononic  heat  transfer  not  experienced  in  the  shorter  ones. 

The  effect  of  chirality  on  A  is  evaluated  by  testing  CNTs  with 
D  ^4.5  nm  and  different  chirality  ((33,33)  armchair;  (43,23)  chiral 
and  (57,0)  zigzag)  as  well  as  D  ^  5.0  nm  ((37,37)  armchair;  (47,27) 
chiral  and  (64,0)  zigzag),  being  the  length  fixed  at  20  nm.  Results  in 
Fig.  4c  show  that  armchair  and  chiral  CNTs  have  similar  A  (average 
difference  below  1%),  whereas  zigzag  CNTs  have  a  slightly  lower  A 
than  armchair  and  chiral  ones  (about  15%  less).  Note  that  the  A  ~  D 
trend  is  confirmed.  While  a  similar  proportionality  is  also  found  in 
other  works  [79],  in  literature  there  is  no  general  consensus  on 
effect  of  chirality  on  A.  However,  phonon  density  of  states  is  almost 
not  influenced  by  CNT  chirality  [80],  which  also  means  that 
chirality  has  a  negligible  effect  on  A. 

Thermal  conductivities  of  SWCNTs  are  compared  to  thermal 
conductivities  of  CNTs  with  same  length,  chirality  and  diameter  but 
being  part  of  DWCNTs,  as  inner  or  outer  tubes.  The  analyzed  cases  are 
(from  left  to  right  in  Fig.  4d):  10  nm  long  (5,5)  SWCNT  vs.  inner  tube 
of  (5,5)(10,10)  DWCNT;  10  nm  long  (10,10)  SWCNT  vs.  outer  tube  of 
(5,5)(10,10)  DWCNT;  20  nm  long  (37,37)  SWCNT  vs.  outer  tube 
of  (33,33)(37,37)  DWCNT;  20  nm  long  (47,27)  SWCNT  vs.  outer  tube 
of  (43,23)(47,27)  DWCNT;  20  nm  long  (64,0)  SWCNT  vs.  outer  tube  of 
(57,0)(64,0)  DWCNT.  As  evident  in  Fig.  4d,  a  general  slight  increase  in 
A  (about  20%  more  in  average)  is  noticeable  when  the  CNT  is  part  of  a 
DWCNT,  regardless  of  the  fact  that  it  is  the  inner  or  the  outer  tube.  To 
our  knowledge,  this  effect  has  never  been  reported  before.  However, 
the  observed  trend  may  be  due  to  the  stress  induced  on  the  inner/ 
outer  DWCNT  tubes  by  reciprocal  Van  der  Waals  interactions,  which 
has  been  demonstrated  to  enhance  A  [81  . 

Once  geometrical  effects  on  A  are  explored,  other  simulations 
have  been  performed  for  estimating  the  thermal  conductivity  of 
chemically  functionalized  SWCNTs.  In  this  case,  parallel  overlap¬ 
ping  CNTs,  covalently  bonded  each  other  by  a  variable  number  of 
C-O-C  joints,  are  considered  (see  Fig.  le).  To  this  purpose,  as 
shown  in  Fig.  lc,  carbon  rings  are  thermostated  (green  dotted  area) 
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Fig.  4.  Effect  of  CNTs  geometry  on  thermal  conductivity:  (a)  effect  of  SWCNT  diameter;  (b)  effect  of  SWCNT  length;  (c)  effect  of  SWCNT  chirality;  and  (d)  comparison  between 
thermal  conductivities  of  CNTs  in  single-walled  or  double-walled  configurations. 
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at  the  ends  of  the  lower  (5,5)  SWCNT  (L=  10  nm,  with  an  overlap 
between  the  two  jointed  SWCNT  of  2.5  nm  along  the  CNT  axis),  in 
order  to  extract  temperature  profile  and  heat  flux  between  the 
thermostats.  Cases  with  0,  1  and  4  C-O-C  joints  are  considered, 
and  the  resulting  X  are  76,  60  and  60  Wm-1  K-1,  respectively. 
Again,  X  is  75  W  m-1  K-1  in  the  sole  SWCNT  case,  as  previously 
shown.  Therefore,  although  X  of  SWCNTs  is  not  affected  by  overlap 
(case  0  joints);  a  significant  drop  in  X  (about  20%)  is  found  by 
adding  one  surface  functionalization  to  SWCNT  (case  1  joint). 
Moreover,  the  conductivity  X  undergoes  the  same  decrease  by 
slightly  varying  the  amount  of  functionalizations  (case  4  joints). 

As  a  matter  of  fact,  Oxygen  joints  induce  rehybridization  of  the 
functionalized  carbon  atoms,  whose  structure  changes  from  sp2  to 
sp3.  Hence,  the  presence  of  sp3  carbons  in  the  crystalline  structure 
of  CNT  causes  a  localized  phonon  scattering,  which  reduces  the 
phonon  mean  free  path  thus  lowering  X  up  to  70%,  according  to  the 
density  of  functionalizations  [82  .  Note  that,  the  arrangement  of 
defects  has  also  a  great  importance  in  controlling  X:  Once  the  first 
defect  is  introduced,  patterned  functionalizations  with  stripes 
parallel  to  heat  flux  (as  our  case)  lower  X  more  gradually  than 
perpendicular  or  random  ones,  while  the  amount  of  functionaliza¬ 
tions  is  increased  [83]. 


4.2.  Thermal  boundary  resistance  in  carbon  nanotube  networks 


The  localized  sources  of  thermal  boundary  resistance  are 
identified  in  the  simulated  setups,  which  mimic  the  building  blocks 
of  more  complex  CNT  network.  Rk  within  and  between  DWCNTs  is 
then  evaluated,  and  the  effect  of  CNT  overlap  and  bonding 
compared. 

In  general,  the  thermal  boundary  resistance  due  to  local  phonon 
scattering  is  calculated  as 


where  AT  is  the  time  averaged  temperature  discontinuity  at  the 
considered  interface.  Let  us  consider  the  geometry  in  Fig.  2d  as 
representative  of  the  building  block  of  CNT  networks:  Thermal 
boundary  resistances  can  be  found  between  the  layers  of  DWCNTs 
(ft/c i  and  Rk2 ,  also  known  as  intralayer  Kapitza  resistances)  as  well 
as  between  overlapping  adjacent  CNTs  (ftk> 12,  also  known  as 
interlayer  Kapitza  resistance).  In  the  first  case  AT  refers  to  the 
time  and  space  averaged  temperature  drop  at  steady  state  between 
inner  and  outer  CNTs;  while  in  the  second  case  AT  is  the  time  and 
space  averaged  temperature  drop  at  steady  state  between  over¬ 
lapping  area  of  parallel  CNTs.  In  both  cases,  cp  is  induced  by  the 
different  temperatures  of  thermostats  at  the  ends  of  the  inner  CNTs 
(Fig.  lc  and  d,  red  highlighted  carbon  rings).  Note  that  a  local 
temperature  drop  is  also  observed  at  the  interface  between 
thermostated  and  remaining  carbon  atoms;  however,  the  corre¬ 
sponding  Rk  (about  2  x  10-11  m2KW-1)  can  be  2-3  orders  of 
magnitude  lower  than  inter-  or  intralayer  thermal  boundary 
resistances  and  thus  neglected  [84  . 

First,  ft/u  and  Rk2  are  evaluated  within  a  pair  of  parallel  (5,5) 
(10,10)  10  nm  long  DWCNTs  with  a  2.5  nm  overlap  and  a  number  of 
joints  ranging  from  0  to  7  (Fig.  2d).  While  no  significant  change 
in  the  value  of  ft,u  or  Rk2  is  noticed  by  varying  the  amount 
of  O-joints,  Rk2  (i.e.  the  Kapitza  resistance  within  DWCNT  with 
280  K  thermostat  attached  to  the  inner  tube)  is  found  to  be  slightly 
higher  than  Rk  l  (i.e.  Kapitza  resistance  within  DWCNT  with  320  K 
thermostat  attached  to  the  inner  tube),  namely  (2.6  +  0.3  s.d.)x 
10” 10  m2  KW'1  and  (2.4  ±  0.4  s.d.)x  10”10m2KW”\  respec- 
tively.  This  difference  is  due  to  the  inverse  proportionality  of  Rk 
with  temperature  [85  . 

Second,  the  interlayer  Kapitza  resistance  is  studied  by  changing 
the  overlap  between  SWCNTs  and  the  amount  of  oxygen  joints 
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Fig.  5.  Effect  of  the  overlap  a  on  the  interlayer  Kapitza  resistance  Rk  u  for  a  couple 
of  parallel  (5,5)  SWCNTs. 


between  adjacent  SWCNTs  or  DWCNTs.  The  overlap  effect  on  Rk  u  is 
investigated  by  varying  the  length  of  the  overlapping  region  between 
two  parallel  (5,5)  SWCNTs  (L=10  nm  both  tubes;  ft  =  0.23  nm;  a= 1, 
2.5,  5,  7.5  and  10  nm),  as  sketched  in  Fig.  5.  Rk  u  is  found  to  follow  an 
exponential  trend  as  a  function  of  the  overlap  length  (Fig.  5),  ranging 
from  8.0  x  10”9  m2I<W”1  (a=lnm)  to  1.3  x  10”9  m2  K  W”1 
(a  =  10  nm),  i.e.  a  6-fold  reduction,  which  can  be  fitted  by 

Rk  =  'Ecie-^x+R°.  (6) 

i 

The  asymptotic  interlayer  Kapitza  resistance  Rk  is  found  to  be 
1.4  x  10 ~9  m2  KW"1,  and  MD  results  are  best  fitted  (ft2  >  0.99) 
by  a  single  exponential  decay  (Eq.  (6)  with  i=l,  x  =  a,  C\  =  13. 9x 
10-9  m2  K  W-1,  /T,  =0.7  nm-1).  The  remarkable  decrease  in  Rk  u 
is  due  to  the  increase  of  the  available  surface  for  heat  transfer 
between  the  contiguous  CNTs  through  van  der  Waals  interactions. 
Even  for  different  setups,  similar  reductions  of  Rku  with  the 
increasing  overlapping  area  between  CNTs  have  been  reported 
[84,86-88].  Moreover,  previous  experimental  studies  agree  with 
the  obtained  Rku  values:  For  instance,  Yang  et  al.  report  experi¬ 
mental  thermal  boundary  resistances  between  MWCNTs  around 
1.0  x  10-9m2I<W-  [89  .  Finally,  the  effect  on  Rku  of  adding 

covalent  bonds  between  a  couple  of  parallel  (5,5)  SWCNTs 
(L=10  nm  in  both  tubes;  h  =  0.23  nm;  a  =  2.5  nm;  number  of  join- 
ts  =  0,  1,  4,  8)  or  (5,5)(10,10)  DWCNTs  (L=10nm  in  both  tubes; 
/i  =  0.23nm;  a=2.5  nm;  number  of  joints  =  0,  1,  4,  7)  is  also 
investigated.  In  the  SWCNTs  setup  in  Fig.  6,  AT0j  is  the  average 
temperature  difference  between  the  top  and  bottom  tube  in  the 
overlapping  area  of  the  configuration  with  0  joints;  whereas  AT4j 
refers  to  the  4  joints  case.  In  both  cases,  the  interlayer  Kapitza 
resistance  Rk  u  is  evaluated  by  means  of  Eq.  (5)  and  using  AT0j  and 
AT4j,  respectively.  In  the  DWCNTs  setup  in  Fig.  7,  instead,  bigger 
dots  represent  inner  tubes  (thermostats  are  attached  at  the  tube 
ends)  while  smaller  diamonds  are  representative  of  the  outer 
tubes;  as  an  example,  A Tin_out  is  the  average  temperature  differ¬ 
ence  between  inner  and  outer  tubes  of  top  DWCNT  in  the  0  joints 
case,  which  allows  to  evaluate  the  intralayer  Kapitza  resistance  Rk  l 
by  means  of  Eq.  (5).  As  shown  in  Figs.  6  and  7,  the  introduction 
of  C-O-C  joints  causes  a  decrease  of  the  average  temperature 
difference  between  contiguous  SWCNTs  or  DWCNTs  (e.g. 
AT4j  <  Ar0j  if  the  SWCNTs  case  with  4  joints  is  compared  to  the 
0  joints  one,  as  depicted  in  Fig.  6)  thus  reducing  significantly  the 
value  of  Rku.  If  the  SWCNTs  setup  (as  shown  in  Fig.  2c)  is 
considered,  Rku  ranges  from  3.6  x  10-9  m2I<W-  (0  joints)  to 
0.7  x  10"9m2  KW-1  (8  joints),  i.e.  a  5-fold  reduction;  whereas  in 
the  MWCNTs  setup  (as  shown  in  Fig.  2d)  Rku  ranges  from 
1.9  x  10"9  m2  KW"1  (0  joints)  to  0.4  x  10-9  m2  KW"1  (8  joints), 
i.e.  a  5-fold  reduction.  In  both  cases,  Rku  shows  an  exponential 
decay  by  increasing  the  amount  of  joints  (Fig.  8),  which  is  again 
fitted  with  ft2  >  0.99  by  Eq.  (6),  by  considering  i= 2  (double 
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decaying  exponential  function  with  C\  =2.0x10  9m2I<W  \ 
=4.7,  c2  =  1.0  x  10~9  m2  KW'1  and /32=0.3  in  case  of  SWCNTs; 
ci  =1.1  x  10~9  m2  KW"1,  /?a  =2.3,  c2  =  0.7  x  10~9  m2  KW"1 

and  /?2  =  0.1  in  case  of  DWCNTs)  and  x  as  the  number  of  joints. 
The  asymptotic  interlayer  Kapitza  resistance  Rk  (i.e.  the  minimum 


Fig.  6.  Steady  state  temperature  profiles  of  a  pair  of  overlapping  SWCNTs 
connected  by  0  (black  dots)  or  4  (yellow  dots)  C-O-C  covalent  bonds  (setup  shown 
in  Fig.  lc).  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the 
reader  is  referred  to  the  web  version  of  this  article.) 
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Fig.  7.  Steady  state  temperature  profile  of  a  pair  of  overlapping  DWCNTs  connected 
by  0  (black  dots)  or  4  (yellow  dots)  C-O-C  covalent  bonds  (setup  shown  in  Fig.  Id). 
(For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is 
referred  to  the  web  version  of  this  article.) 


V  -3*  M-n  •  +  ^ 

f 

-  ^  W\ 

...  i.  .**  •• 

^fin-out 

♦**  *♦♦***♦♦♦*♦*%♦♦ ~  j 

i  * 

i 

i 

i 

i 

i 

i 

■  •  0  joints  (In) 

1 

1 

i 

i 

i 

i 

•  0  joints  (Out) 

4  joints  (In)  j 

'  4  joints  (Out) 

"“'Initial  state 

111,1 

i 

Fig.  8.  Interlayer  Kapitza  resistance  Rk  u  as  a  function  of  the  number  of  C-O-C 
covalent  bonds  (i.e.  joints)  between  a  couple  of  parallel  (5,5)  SWCNTs  or  (5,5) 
(10,10)  DWCNTs. 


value  of  the  Kapitza  resistance  achievable  by  adding  a  large  amount  of 
C-O-C  joints  between  the  CNTs)  is  found  to  be  0.5  xlO-9  m2  KW' 
for  SWCNTs,  and  0.1  x  10“ 9  m2  K  W_1  for  MWCNTs.  The  reported 
reduction  in  thermal  boundary  resistance  is  due  to  the  supplementary 
bonded  pathways  -  in  addition  to  non-bonded  ones  -  for  carrying 
thermal  vibrations  (i.e.  phonons)  between  adjacent  CNTs.  DWCNTs  and 
SWCNTs  show  the  same  decaying  exponential  trend  of  Rku,  but 
DWCNTs  values  are  shifted  down  because  of  their  larger  contact  area 
in  this  case  (i.e.  lower  Rku  due  to  overlap  effect),  which  gives  rise  to 
the  observed  offset.  As  in  the  above  analysis  on  the  overlap  length,  the 
decrease  in  Rk  u  due  to  the  increase  of  the  number  of  links  between 
CNTs  is  in  general  agreement  with  both  experimental  [90,91]  and 
modeling  studies,  even  with  different  MD  approaches  and  setups 
[57-60,92]. 


5.  Discussion 

According  to  the  published  results  and  the  performed  MD 
simulations,  some  key  factors  controlling  thermal  transmittance 
of  CNT  networks  are  highlighted. 

First,  CNTs  should  be  a  few  hundreds  of  nanometer  long,  for 
allowing  the  existence  of  phonons  with  longer  wavelengths  thus 
increasing  A  but  avoiding  technological  problems  such  as  structure 
stability  and  presence  of  impurities,  which  are  typical  for  micro¬ 
meter  long  CNTs. 

Second,  an  ordered,  parallel  and  partially  overlapping  arrange¬ 
ment  of  CNTs  should  be  pursued  within  the  network,  though 
avoiding  direct  contact  between  CNTs,  which  causes  defects  in 
their  crystalline  structure  and  then  a  decrease  of  their  internal  A 
[93].  Note  that,  because  of  the  exponential  decay  governing  the 
relation  between  Rk  and  the  overlap  length,  even  a  short  overlap  is 
highly  beneficial  to  the  overall  thermal  transmittance. 

Third,  the  introduction  of  a  few  short  C-O-C  joints  between  the 
CNTs  constituting  the  carbon-based  network  significantly  improves 
the  thermal  transmittance  of  the  composite.  As  also  noticed  in 
other  studies  [59,60],  short  and  stiff  bonds  (such  as,  C-O-C  or 
C-CH2-C)  chemically  affine  to  C-C  should  be  preferred  for  reducing 
Rk,  because  of  the  corresponding  shorter  distance  between  CNTs, 
which  allows  a  higher  heat  transfer  via  Van  der  Waals  interactions, 
and  the  coupling  of  vibrational  modes  of  the  nanoconstructs 
within  the  network. 

Therefore  we  would  like  to  stress  that,  because  of  the  decaying 
exponential  trend  of  Rk  with  respect  to  the  number  of  joints,  even  a 
single  bonded  interaction  could  be  sufficient  for  a  significant 
improvement  of  thermal  transmittance. 

Finally,  the  problem  of  finding  a  balanced  number  of  bonded 
joints  between  CNTs  is  also  addressed,  because  the  functionaliza¬ 
tion  of  nanotubes  disturbs  pure  sp2  bonds  thus  inducing  lower 
thermal  conductivity.  Therefore,  when  high  thermal  performances 
of  CNT  networks  are  demanded,  the  best  design  compromise 
between  the  introduction  of  functionalizations  (i.e.  phonon  scat¬ 
tering  sources)  and  the  need  for  reducing  the  thermal  boundary 
resistance  is  to  add  a  limited  amount  of  short  joints  between  CNTs, 
chemically  affine  to  C-C  bonds.  Hence,  the  quality  and  quantity  of 
CNTs  interlinks  as  well  as  the  aligned  and  partially  overlapping 
arrangement  of  CNTs  appear  to  be  crucial  for  improving  thermal 
conductance  of  CNT  networks  and  their  composite  materials,  even 
at  low  filler  content. 


6.  Conclusions 

In  this  work,  we  have  reviewed  and  provided  numerical  evidence 
of  the  impact  of  geometrical  parameters  and  chemical  functionali¬ 
zations  on  the  overall  thermal  transmittance  in  large  complex 
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networks  made  up  of  carbon  nanotubes  (CNTs).  The  considered 
geometrical  parameters  are  mainly  dictated  by  the  nanotube  size, 
while  functionalizations  were  implemented  by  the  introduction  of 
short  C-O-C  interlinks  between  two  adjacent  nanotubes. 

For  the  sake  of  simplicity,  we  have  focused  on  the  building 
blocks  of  the  above  networks,  so  that  the  considered  configurations 
consist  of  single  and  pairs  of  both  single-walled  carbon  nanotubes 
(SWCNTs)  and  double-walled  carbon  nanotubes  (DWCNTs)  with 
different  dimensions  and  number  of  interlinks. 

We  observed  that  as  far  as  single  tubes  are  concerned,  diameter 
has  a  weak  influence  on  the  thermal  conductivity  A,  whereas  a 
more  significant  impact  is  found  instead  for  the  length.  Apart  from 
secondary  effects,  the  above  conclusion  also  holds  for  each  nano¬ 
tube  within  a  nanotube  aggregate.  In  the  latter  case,  high  thermal 
boundary  resistances  (interlayer  Kapitza  resistance)  are  observed 
for  contiguous  nanotubes  since  the  interaction  between  particles 
only  relies  upon  secondary  forces  (van  der  Waals),  and  it  exponen¬ 
tially  decays  to  an  asymptotic  value  by  increasing  the  overlapping 
region.  Moreover,  a  further  reduction  in  the  interlayer  Kapitza 
resistance  can  be  observed  owing  to  the  presence  of  interlinks 
between  two  nanotubes.  Here,  in  particular,  a  double  exponential 
decay  proves  to  best  fit  the  measured  values.  Our  results  show  that 
few  connections  (two  or  three)  are  already  sufficient  to  produce  a 
significant  improvement,  while  a  further  increase  in  the  number  of 
joints  only  produces  a  marginal  reduction  of  the  above  resistance. 

The  present  work  may  provide  useful  guidelines  in  designing 
the  process  for  fabricating  novel  composite  materials  of  great 
interest  in  modern  thermal  sciences,  namely  thermal  storage  and 
polymeric  heat  exchangers. 
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